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Abstract 

We consider an eigenvalue problem for a divergence form elliptic operator with high contrast pe- 
riodic coefficients with period e in each coordinate, where e is a small parameter. The coefficients are 
perturbed on a bounded domain of 'order one' size. The local perturbation of coefficients for such operator 
could result in emergence of localized waves - eigenfunctions with corresponding eigenvalues lying in the 
gaps of the Floquet-Bloch spectrum. We prove that, for the so-called double porosity type scaling, the 
eigenfunctions decay exponentially at infinity, uniformly in e. Then, using the tools of two-scale conver- 
gence for high contrast homogenization, we prove the strong two-scale compactness of the eigenfunctions 
of Ag. This implies that the eigenfunctions converge in the sense of the strong two-scale convergence to 
the eigenfunctions of a two-scale limit homogenized operator Aq, consequently establishing 'asymptotic 
one-to-one correspondence' between the eigenvalues and the eigenfunctions of these two operators. We also 
prove by direct means the stability of the essential spectrum of the homogenized operator with respect to 
the local perturbation of its coefficients. That allows us to establish not only the strong two-scale resolvent 
convergence of A^ to Aq but also the Hausdorff convergence of the spectra of A^ to the spectrum of Aq, 
preserving the multiplicity of the isolated eigenvalues. 

Keywords: localized modes, elliptic operators, perturbed periodic operators, multiscale methods, two- 
scale convergence, high-contrast homogenization 

AMS Subject Classifications: 35B27, 35P99 

1 Introduction 

In this paper we consider a high contrast two-phase periodic medium with a small period and with a 'finite 
size' defect filled by a third phase, see Fig. [H This physically represents, for instance, a simplified model 
of cross-section of a photonic crystal fiber. Mathematically, the problem relates to a compact perturbation 
of e-periodic coefficients in a divergence form elliptic operator A^. The behaviour of and its spectral 
characteristics as e — > is of the main interest. A similar problem is considered in p^J using the method of 
asymptotic expansions, but the present study pursues different aims and approaches the problem from another 
direction, namely developing an appropriate version of the two-scale convergence technique [T9l|2l|2Tj. As 
a result we obtain a complete description of the asymptotic (with respect to e) behaviour of the localized 
modes and other spectral characteristics for the operator Af, in terms of an explicitly described (two-scale) 

*The author thanks V.P. Smyshlyaev and I.V. Kamotski for their help and attention to this work. The author also thanks 
V.V. Kamotski for valuable suggestions. 



limit operator Aq. For other recent applications of the high contrast homogenization techniques see also 

ini El El El [mil]. 

In the absence of a defect, Zhikov considers in [22] a divergence form elliptic operator (denoted by A^ in 
[22]) with periodic coefficients corresponding to a double-porosity model [31 [9] (^^ in our notation is obtained 
from A^ by a compact perturbation of its coefficients). Operators of such type have the Floquet-Bloch essential 
spectrum, displaying a band-gap structure. Zhikov proves that the spectra of A^ converge in the sense of 
Hausdorff to the spectrum of a certain two-scale homogenized operator with constant coefficients, see also 
[14:\ I21|. and that Aq is the limit of Af, in the sense of strong two-scale resolvent convergence. The spectrum 
of ^0 is purely essential and displays an explicit band-gap structure. It is well known, see e.g. [20l[l3], that 
in the case of a compact perturbation of periodic coefficients in the elliptic operator A^ its essential spectrum 
remains unperturbed. The only extra spectrum that can emerge in the gaps due to the perturbation is a 
discrete one (isolated eigenvalues with finite multiplicity) Such an extra spectrum does emerge at least 
under some assumptions, e.g. [13^ I16j. The latter corresponds physically to localized modes emerging near 
the defect. In order to establish the strong two-scale convergence of the eigenfunctions of A^ we need their 
strong two-scale compactness. The latter requires in turn an exponential decay of the eigenfunctions uniform 
in £. 

The problem of wave localization (i.e. of the existence of eigenvalues with corresponding eigenfunctions 
decaying exponentially) in the gaps of the essential spectrum has been intensively investigated for a wide 
range of differential operators over the last decades. The results obtained up to date ensure the exponential 
decay of eigenfunctions of A^ for a fixed e, see e.g. |13j . However this is insufficient for establishing the 
required compactness. Moreover, the developed methods, e.g. [5] and [l3] (the latter using the method of 
Agmon[T]), seem to be insufficient for the present purpose. The reason is that in order to obtain the uniform 
exponential decay one has to perform some kind of two-scale asymptotic analysis, investigating the behaviour 
of the eigenfunctions on small and large scales simultaneously. To achieve this we supplement the method of 
[1] by the related two-scale techniques, which play a crucial role. As a result, we obtain a uniform estimate 
with the decay exponent a (see p.ip and (|2.14p below) which ensures the compactness, but may also be of 
an independent interest. On one hand, it is sharp in a sense. On the other hand, it behaves qualitatively 
entirely different compared to e.g. the one in [5]: while the one in [5] is proportional to the square root of 
the distance to the gap end, the decay exponent we derive becomes large on approaching the left end of the 
gap and small near the right end. 

The structure of the paper is the following. We first define the problem in Section [21 describe the two-scale 
limit operator and state the main result. We then consider a subsequence of eigenvalues of A^ converging 
to some point Aq lying in a gap of the spectrum of Aq. In Section [3] we prove (Theorem 13. ip the uniform 
exponential decay for the eigenfunctions of A^. Section [His devoted to the proof of a main auxiliary lemma 
that is employed in the previous section, which may also be of an independent interest. In Section [5] we list 
some properties of the two-scale convergence and several related statements which we use in the next section. 
Employing the uniform exponential decay, we establish in Section [6] (see Theorem 16. ip the strong two-scale 
compactness of (normalized) eigenfunctions of A^, see e.g. [211 [22]. This implies that, up to a subsequence, 
the eigenfunctions two-scale converge to a function, which is eventually proved to be an eigenfunction of the 
two-scale limit operator with a defect, which could be considered as a perturbation of Aq. Accordingly Aq 
is an eigenvalue of ^o- The two-scale convergence of the eigenfunctions together with the results of [16] on 
the existence of the eigenvalues in the gaps and related error bounds allow us to make a conclusion about the 
'asymptotic one-to-one correspondence' between eigenfunctions and eigenvalues of the operators A^ and 
as e — > 0. In the last section we prove by direct means (via the Weyl's sequences) the stability of the essential 
spectrum of Aq with respect to the local perturbation of its coefficients (see Theorem 17. ip . Thereby this 

^We do not concern in this paper the issue of whether embedded eigenvalues can emerge on the bands as a result of the 
perturbation. 
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Figure 1: A defect in a rapidly oscillating high contrast periodic medium, cf. [16, Fig. 1]. 



establishes the convergence of the spectra of to the spectrum of Aq in the sense of Hausdorff (Theorem 
O). 



2 Notation, problem formulation, limit operator and the main result 

We will use the following notation for the geometric configuration visualized on Figure [H cf. . Consider 
a periodic set of unit cubes 

{Q:Q = [0,ir+C,eeZ"}. (2.1) 

Let Fq be an open periodic set with period one in each coordinate such that Fq Q <^ Q is a connected 
domain with infinitely smooth boundary. We denote FqCiQ hy Qq and Q\Qq by Qi. Notice that the position 
of the particular set Qo, Qi or Q depends on ^ E Z", however we will not reflect this fact in the notation to 
simplify the latter. Regularity assumptions on the boundary could be relaxedH Let ^2 be a bounded domain 
with a sufficiently smooth boundary, containing the origin; its complement is denoted by J^i, = M"\r22. 
We define the 'inclusion phase' or the 'soft phase' Oq as 



where e > is a small parameter. The set of inclusions eQq which intersect the boundary of ^2 is denoted by 

Oq. The 'matrix phase', denoted by J7f , is the complement to the inclusions in f^i, i.e. fif = r2i\(i7Q U J^q). 

'Defect domain' 01 is defined by i72\^o- notation 9q for the characteristic function of a set 

and Bji for the open ball of radius R centered at the origin. 
We consider an elliptic operator A^, self-adjoint in (K"), 



AeU' 



V • ( a{x, e)Vn^(x) ) , x G 



in particular the eigenvalue problem 



(2.2) 



(2.3) 



^In particular, the results on the two-scale convergence stated in the paper remain valid at least under the assumption of 
Lipschitz regular boundaries. The e^'^'^-order bounds, as they were obtained in [16], require higher regularity. 
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for its point spectrum. The coefficient a{x, e) is given by the formula 

a{x, e) = 

where measurable ao{x,£) is such that 





X 




ai, 


X 




02, 


X 




ao(x,e), 


X 




for all 




or ao{x, £ 



(2.4) 



oo for all e. (2.5) 



Here oq, ai, 02, ^o, -^o and 6 are some positive constants independent of e, 6* E (0,2]. Notice that this 
includes as particular cases e.g. the case of 'removed' boundary inclusions, i.e. a{x,e) = ai if x € fig l^i, 
a(x, e) = 02 if X E Oq n $72, and the case of the 'full' inclusions, ao(x, e) = oq e^. The domain of is defined 
in a standard way via Friedrichs extension procedure with a bilinear form, see ^Mi below, defined on H^{W). 

For any e > the operator Af, is an operator with e-periodic coefficients, which are compactly perturbed 
(within bounded domain f^l U fip). This implies (e.g. [20^ I13j) that its essential spectrum coincides with the 
Floquet-Bloch spectrum of the associated 'unperturbed' operator A^, with only extra spectrum being hence 
the discrete spectrum in the gaps of lei Note that the spectrum of Af. contains gaps for small enough e, cf. 
[m [2TI [22], and there is often an extra discrete spectrum in the gaps of aessi-^e)-, e.g. [E]. By definition, 

E -ff^(M"), ^ 0, is an eigenfunction of the eigenvalue problem ()2.3p with an eigenvalue if 

j a{x,£)Vu^ -Vw dx = \e j u'^wdx (2.6) 

for all w E H\R''). 

The aim of this work is to establish that as e — s- the operator A^ converges in the appropriate sense 
(namely, in the sense of two-scale convergence, see Section [5]) to a 'two-scale' limit operator Aq, which we 
describe next. For the rest of the present section we assume that Q = [0, 1)", considering all functions of 
two variables (x, y) to be 1-periodic in each coordinate with respect to y. The 'two-scale' limit operator ^0 
is analogous to the one introduced in the defect free setting by Zhikov [21^ [22] and acts in a Hilbert space 



Ho ■■= \u{x,y) E (M" x Q) 



u{x, y) = uq{x) + v{x, y),uo e 



L^W'),ve L^nr,L\Qo))'^, (2.7) 



with the natural inner product inherited from L^(M" x Q) and TIq being its closed subspace, cf. [22j. It 
is implied that v is extended by zero for y Qi or x & i^2- The operator ^0 is defined as generated by a 
(closed) symmetric and bounded from below bilinear form Bq{u,w) acting in a dense subspace 

V = (M") + {ni,H^{Qo)) (2.8) 

of Ho = L"^ (M") -I- (ill, L^(Qo)) , which is defined as follows: ioi u = uq + v ,w = wo + z ^ V, 

Bo{u, w) = a2 j Vuo • Vwo dx + j A^"^Vuo ■ Vwq dx + qq j j VyV ■ VyZ dy dx. (2.9) 
Here A^""^ = (Af""^) is the standard ' 'porous" homogenized (symmetric, positive-definite) matrix for the 



periodic medium as described above but when no defect is present and with gq = 0, see e.g. [151 §3-1]: 

= inf / ail^ + Vw\^ dy (^ E M") . (2.10) 



This does not rule out possible emergence of embedded eigenvalues on the bands, not considered in this paper. 
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Figure 2: /3(A), cf. 



Here C^^{Q) stands for the set of infinitely smooth functions with periodic boundary conditions. Then one 
can see (cf. |22j ) that the form is indeed bounded from below, densely defined and closed. Hence, according 
to the standard Friedrichs extension procedure, e.g. |20j, can be defined as a self-adjoint operator with a 
domain ^^{Aq) C V. A function u^{x,y) = uo{x) + v{x,y) G V, u^{x,y) ^ 0, is an eigenfunction of the limit 
operator Aq corresponding to an eigenvalue Aq if and only if 



B(){u^,w) = Xo / I {uq + v){wq + z) dy dx. 



(2.11) 



for any w = wq + z ^ V (we assume where it is possible that a function defined on a smaller domain is 
extended by zero on a larger domain) 

The 'unperturbed' operators Agr and Aq could be defined analogously to Af, and Aq formally setting above 
1^2 = and Q.I = M". (See also [2T1[22], where these operators are denoted by yl^ and A respectively.) 

We next describe a function /3(A) which was introduced by Zhikov [21^ [22] (cf. also [8j) and plays 
an important role in our considerations. Let \j and (^j, j = 1,2,..., be eigenvalues and corresponding 
orthonormalized eigenfunctions of operator T defined as 



Tf := -aoA/, / G hHQq) n H\Qq). 



(2.12) 



Note that the eigenvalues of T belong to the spectrum of Aq, see [21]. For X Xj, j > 1, denote by b the 
solution to 

Tb-Xb = -aoAb -Xb = l, 6 G H^{Qo). (2.13) 



The function /3(A) is defined by 



/3(A) := X{1 + X{b)y) = X + X^J2 



■ 1 - ^ 



(2.14) 



where {f)y := / f{y) dy. It is well-defined for any A except A = Aj with {ipj)y / 0, monotonically increasing 

Q 

between such points, see Figure [2J This function describes the structure of a{AQ), see [21]. Namely, the 
intervals where /3(A) > correspond to the bands of the spectrum of Aq. Isolated points of the spectrum of 



^Explicit example in [161 §5] ensure the existence of isolated eigenvalues of Aq of finite multiplicity in the gaps of g{Ao) in a 
particular situation. 
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i-e. Aj such that {'fj)y = and /3(Aj) < 0, can also be regarded as degenerate bands. The intervals on 
which /3(A) < (excluding Aj) are gaps. 

It was shown in [22] (see also 121] ) that (t{Ai,) converges in the sense of Hausdorff to a{AQ), while 
Af, converges to Aq in the sense of the strong two-scale resolvent convergence (cf. Sections [5] and [6] below) 
implying the convergence of spectral projectors, etc. 

We aim at showing that similar as well as some further results hold for the perturbed operators. Namely, 
our main result is the following 

Theorem 2.1. The operator A^ converges to Aq in the sense of the strong two-scale resolvent convergence. 
Hence the spectral projectors also strongly two-scale converge away from the point spectrum of Aq. The 
spectrum of A^ converges in the sense of Hausdorff to the spectrum of Aq. Let Aq be an isolated eigenvalue of 
multiplicity m of the operator Aq in the gap of its essential spectrum. Then, for small enough e, there exist 
exactly m eigenvalues X^^i of A^ (counted with their multiplicities) such that 

|Ae,i-Ao| <Cei/2^i = l,...,m, (2.15) 

with a constant C independent of eJl If for some sequence Sk a sequence of eigenvalues A^^, of A^ 
converges to Xq which is in the gap of the essential spectrum of Aq, then, Xq is an isolated eigenvalue of Aq 
of a finite multiplicity m and for large enough k, Ag^. G {Aej.,i, i = 1, . . . , rn}. 

A key part in establishing the latter is in controlling the behaviour at infinity of the eigenfunctions 
corresponding to the extra point spectrum which may appear in the spectral gaps of the unperturbed operator. 
A central property providing this is a uniform exponential decay of the eigenfunctions which we prove next. 

3 Uniform exponential decay of the eigenfunctions of 

Let Aq be a point in a gap of a{AQ), i.e. such that /3(Ao) < and Aq ^ Xj for all j. Assume Aq is an 
accumulation point of the point spectra of A^, i.e. for some subsequence — > there exist eigenvalues A^^. 
of Ag such that A^^, — > Aq as /c — > oo. (Notice that the results of [I3l[l6] ensure in particular that such series 
do exist.) We formulate the main result of this section (and also one of the principal results of the paper) in 
the following statement. 

Theorem 3.1. Let X^^. and n^'-' be sequences of eigenvalues of the operator A^ and corresponding eigen- 
functions normalized in L^(M"), where is some positive sequence converging to zero as /c —> oo. Let Xq 
be such that (3{Xq) is negative and Xq is not an eigenvalue of the operator T given by /12.12\} . Suppose that 
Aej. converges to Xq. Then for small enough eigenfunctions u^^ decay uniformly exponentially at infinity, 
namely, for 

< a < V-/3(Ao)/«i (3.1) 

the following holds: 

||e-l-lu^'=||^2(i,„) <C, 
uniformly in e^, i.e. for any < < ^(q;), with C = C{a) independent of e. 

Proof. We drop the index k in et for the sake of simplification of notation. So, when we say, for instance, 
'sequence A^' we actually mean 'subsequence Ae^.'. 

The plan of the proof is the following. We first derive 'elementary' a priori estimates for the eigenfunction 
outside the set of inclusions JIq U f^Q. Next we study the structure of the eigenfunction at the small scale 

^The error bound (|2.15|l employs the results of [16| requiring, as stated, higher regularity of dQo- The rest of the statement 
of the theorem applies potentially to less regular boundaries. 
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and deduce some vital inequalities for eVu^ inside the inclusions. As a central technical step, we then employ 
in the integral identity (j2.6p a test function with exponentially growing weig ht see (fCTl -dSJSD below, 

and perform some delicate uniform estimates to achieve the result. The main auxiliary technical results are 
proven in Lemma 13.21 and Proposition [TJ 
Step 1. Setting w = in (12. 6p we have 

e^ao\\Vu'\\l2(^^.^ + ai\\Vu^\\l2^^e) + a2\\Vu^\\l2(^^e^^ + \\ao^ix,e)Vu^\\l^^~,^ = Xe\\u'\\l2(^^n) = Xe- 

Therefore 

ll^^llHi(R"\(f2§Uf^o)) ~ ^ ^^''^^ 
uniformly in e. From now on C denotes a generic constant whose precise value is insignificant and can change 
from line to line. 

Step 2. Let us consider the function in a cell eQ corresponding to such ^ = (^(e) G Z", see 12. 11 that the 
corresponding 'inclusion' eQq has a nonempty intersection with Qi. There exists an extension of u'^\eQi to 
the whole cell eQ such that 

Wu^'h^icQo) < C'll^^llL2(eQi), l|V5^llL2(eQo) ^ C\\Vu''\\l2(^^q^), (3.3) 

where C does not depend on e or ^, see e. g. \18\ Ch. 3, §4, Th. 1], which is a version of the so-called 
'extension lemma', see also e.g. [15^ §3.1, L. 3.2]. In particular, we can choose the following extension: 

-V ■ {a{x,e)Vu^{x)) = 0, xen^uhf^, 

which minimizes \\a^^'^{x,e)Vu'^\\]^2(^^Q^'j subject to the prescribed boundary conditions, with ()2.4p and (j2.5p 
ensuring that (13. 3p still holds. From (j3.2p and (|3.3p we conclude that 

<C. (3.4) 

We represent in the form 

u'^{x) = u^{x) + v^{x) (3.5) 

and consider the function G i?g(r2QUJ7Q)|f| In each inclusion sQo C JIqUJIq we have the following boundary 
value problem for v^{x): 

-V • (a(x,e)Vt;=) - A^u^ = AeU^ x G eQo; v^x) =0, x e d{eQo). (3.6) 

When a{x,e) = aoe^, i.e. everywhere in and also in 0.^ in the case ao(x,e) = ooe^, after changing the 
variables x— >y = x/ewe obtain 

- aoAyV^(ey) - XeV^{ey) = XeU%ey), y G Qo, v^{ey) = 0, y G (9(5o- (3.7) 

Since Aq 7^ Aj by the assumptions of the theorem, Ag is separated uniformly from the spectrum of operator 
(j2.12p for small enough e. Hence the resolvent at A^ is bounded uniformly in e and (13. 7p implies 

\\v%ey)\\H^Q,)<C\\u%ey)\\L2(^Q^^y (3.8) 

^In a sense, p.5|) decomposes u'^ into a slowly varying part and rapidly varying . The two are coupled and subsequently 
analyzed simultaneously, which is the essence of two-scale asymptotic analysis. 
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In the case when Aoe"^ ^ < ao{x,e) < Bqe^ ^, G (0,2], we multiply equation ()3.6I1 by and integrate 
by parts to obtain after rescaling 



ao{ey,e)\VyV^ey)\'^dx - Xe / {v%ey)Y dx = / u'' {ey)v^ {ey) dx . (3.9) 



Notice that e aQ{ey,e) > AqE ^ oo as e — > 0. Then using Poincare inequality one easily derives 

E~'\\al/'\7yV^ey)\\l,^Q^) + ||^^(ey)||i.(Q„) < C\\u^iey)\\l.^Qor (3-10) 

for small enough e. Returning in (jS.Sp and ()3.10|) to the variable x we arrive at the following inequality that 
describes the behaviour of v'^ and its gradient in QqUQq, 

||aV2Vt;-(x)||i.(,Q^,) + \\v%x)\\l,^^Q^) < C||5^(x)||i.(,Q^), (3.11) 

with an e- independent constant C. 

Step 3. In order to get the uniform exponential decay of the eigenfunctions we next substitute in ()2.6p 
a test function of a special form: 

w = g^ilxDu^ix). (3.12) 

Here we define function g as follows 

r e"*, te[0,R], 

where R is some arbitrary positive number. The exponent a will be chosen later. This method was employed 
e.g. by Agmon, see [1], but in the present case its realization is not straightforward. Namely, to obtain the 
desired estimates we have to implement the approach of [Ij in the context of the two-scale analysis. We will 
show that g{\x\)u'^ (x), and consequently g{\x\)u^{x), are bounded in L^(M") uniformly with respect to R and 
e. Then we will show via passing to the limit as i? — > oo that we can replace (7(1x1) by e"'^'. 

Remark 1. We cannot use e^"l^lu^(j;) as a test function directly, since it is not known at this stage that this 
function is square integrable. 

The following identity holds by direct inspection 

Vu'Vig'u') = \V{gu')\^ - \Vg\\u'f. (3.14) 
Notice that the absolute value of is bounded by g with a (uniformly in R): 

|V<7(|x|)| <a5(l^l)- (3.15) 
After the substitution of (f3T2]) into (^Ml we have, via ([33]) and (IXTill . 

e'^ao j Vu" ■V{g'^u'')dx + J aoVv^ ■V{g^u^)dx + j a{x,e)\V{gvf)\^ dx- 
-ai j \Vg\^{u^f dx- Xe j g'^{u^fdx-Xejg'^v'u^dx = 
= Xe j g'^u^u" dx + Xe j g'^iu^fdx + j a{x,e)\Vg\'^(u^f dx. 



(3.16) 
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Notice that the right hand side is bounded by some constant C independent of e and R due to ()3.2p . ()3.4p . 
(j3.1ip and the boundedness of the domains of integration. 

We employ (13. 4p , ()3.1ip and the boundedness of Sq to conclude that the second term on the left hand side 
of (j3.16p tends to zero (uniformly in R): 



(3.17) 



as follows. Let us take an arbitrary subsequence vf . Since ||2^||Hi(Kn) is bounded uniformly in e, see ()3.4p . 
the set of functions v!" is weakly compact in H^{Bji), hence strongly compact in L^{Bji) for any R; we take 
R large enough so that CC B^. Then there exists further subsequence 5^ that converges to some function 
"Uo strongly in L'^(Bj^). Then 



lL2(f7f,) 



on < ll^ol 







as Lebesgue measure of the set tends to zero. Since we have chosen in the beginning an arbitrary 
subsequence u^, p.l7p follows. From ()3.1ip and (j3.17p we also obtain 



lL2(Qg) 



0. 



(3.18) 



Step 4. The following Lemma approximates and bounds the last and the first terms (both in a sense of 
a 'two-scale' nature) on the left hand side of (j3.16p . 

Lemma 3.2. There exists Eq > such that for all positive e < £q the following estimates are valid 



Xe J g\'u' dx - (/3(A,) - A,) / g\u'f dx 



< 



(3.19) 



< Ce (||V(5n^)||i2(f,.) + \\gu'\\l2^n.^ni)) + C, 



and 



J Vu'V{g^u') 



dx 



< Ce ( ||V(5n^)||i2(^.) + Ibn^||i2(j,.uf7f) 



■£||2 



+ c 



(3.20) 



where C does not depend on e and R. 



The proof of this lemma is quite technical and we give it in the next section. We make use of Lemma 
and convergence (j3.17p to transform identity (j3.16p into the following inequality, valid for small enough e: 

«i!|V(52")|li2(n^) - ai||(V5)S^||2 2(^,) - /3(A,)||g5^||2 2(f^gu^,)- 
- 25 (||V(5S^)||i.(^.) + ||95^||i2(^guQf)) < C, 

where C is independent of e and R. Notice that [3{Xs) is negative and uniformly bounded away from zero as 
Xe ^ Xq. Applying (j3.15p to the second term on the left hand side we arrive at 



(ai-25)||V(5S^)||i2(o.) + (-/3(A 



a 



^ai -2(5) 115^^11^2 



(f^gUQf) 



(3.21) 
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where 6 > could be chosen arbitrarily small. Hence we should choose a such that — /3(Ao) — a ai is positive, 
i.e. 

a < v^-/3(Ao)/ai. 

Since coincides with e'^'^l on the ball B^, taking 6 small enough and restricting the L^-norms to 

we arrive at 

<C 



uniformly for small enough e. Then passing to the limit as ii ^ oo we obtain 



< C. 



L2(R") 

Step 5. Now we easily get the same estimate for the function n^: 
In each cell we use inequality (13. lip and 



(3.22) 





< 








L2(M") 




L2(M") 



sup e^l^'l < e"v^"e"l^l, Vx G eQ, 

x'eeQ 



to obtain 

and hence, finally, 

uniformly in e. 



LHeQo) 





x\~e 


< C 


^»\x\~e 














< c 






L2(I 







□ 



Remark 2. From (13. 3p . (I3.15P and (I3.2ip it also follows that the gradient of decays exponentially at 
infinity. 



L2(R") 



< c 



(3.23) 



uniformly in e. 

Remark 3. Estimate (13. ip is sharp in a sense. As we will show later, strongly two-scale converges to uq, 



for which y^— /3(Ao)/ai is the optimal estimate for its decay exponent, cf. (|7.15p . 



4 Proof of Lemma [372 • 

Proof. Step 1. First we decompose the function in JIq into the sum of two functions: 

solving the following equations (cf. (|3.7p ): 

- aoAyV^ey) - XeV^ey) = \e{W{ey))y, y G Qo, v^i^y) = 0, y G dQo, 
- aoAyV^{ey) - XeV^{ey) = Xe {vf{ey) - {vf{ey))y) , y G Qo, v^{ey) = 0, y G BQq. 



(4.1) 

(4.2) 
(4.3) 
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The solution of ()4.2p could by presented in the form 

W{ey) = Xe{u%hM. (4.4) 

where is a solution of (j2.13p with A = A^. Due to the uniform (with respect to e) boundedness of the 
resolvent of the operator T in the neighborhood of Aq, the solution of ()4.3p is bounded as follows, 



W{ey)\\m{Q,) < C \\u'{ey) - {u')y\\^^^Q^^ < C W^yu'W^^^^^^^ , 
here we also employed the Poincare inequality. In particular 

\W{x)\\L^ieQo)<^C\\Vu'{x)\\L^eQ), (4-5) 

where C in the inequality does not depend on e or ^ € Z". 

Step 2. At this stage we will need several inequalities which follow from the properties of g and n^. 

Proposition 1. The following estimates are valid for small enough e with constants independent of e and 
the choice of particular eQ: 

¥^"\\meQ) W^^'WlHsQ) < C (||V(55^)||i.(,Q^) + Wgu'Wh^^Q)) , (4.6) 

WWlHsQ) ||V(5'S')|L2(,q) < C (||V(<75^)||i.(,Q^) + \\gu'\\h^,Q)) , (4.7) 

I|V51|l2(,q) ||V(^;V)||^,(^Q) < C (||V(5S^)||i.(,Q^) + \\gu'\\h(^,Q)) . (4.8) 

Proof. Notice that 

supff < e"^^5'(x), X G eQ. 9^ 

eQ ^ ' ' 

We apply ([33]), (IXT51) and (gj]) to get (fO]l : 

= C ll<7n"||i2(,Q) \\^{gu') - (V5)S"||^2(,Q^) < C (||5S"||i2(,Q) ||V(55")||^2(,Q^) + \\gu'\\l2^,Q)) < 

<c(|iV(52^)||i2(,Q,) + ll551li2(,Q)). 

The proof of (j4.7p and (j4.8p is analogous. □ 
Let us show that the entity / g'^v^u'^ dx is relatively small (compared to the first term on the right hand 
side of (j3.19p ). Indeed, applying inequalities (|4.5p and (|4.6p in each cell we obtain 

j g'Wu^dx < ||5'Sli2(,Q„) niL^(.Qo) <Y.'C (i|V(52^)||i2(,Q,) + ||ff51li2(,Q)) . (4.10) 

Considering sets 

[J eQ and (J eQi, 

one can notice that they are "nearly" equal to 

niVJVll and 
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respectively. Namely, 



where 



We introduce two 'correctors' 



and 



QIUQI= [ [j eQ j U fif + \ n 



Of 



eQonfi27^0 

L2(nf _)! 
^1 = Il5^^'^lli2(nf ^uni _)• 



Then inequality (I4.10p transforms into 



Step 3. Now we consider the term J g^v'^u^ dx (cf. p.lOp ) using also (|4.4p and ()2.14p : 

f2S 



A, / g^Wu' dx - (/3(A,) - A,) / g\u'f dx 



2 (~,e\2 



< 



< 



< 



Q 



+ Crl < 



{u')yj {g^u'-{g^u')y)b,dy 
Q 



+ 



{be)y / {g^u'-{g^u')y)u'dy 



+ Ci 



Notice that the mean value of is bounded by its norm in 



I 



dy 



< lk'(^)llL2(Q). 



Similarly, 



{be)y<\\be\\mQo)<C, 
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where C does not depend on e due to the uniform boundedness of (T — A) ^ in the neighborhood of Aq. Via 
the Poincare inequahty we derive 



Q 



< h^u' - {9^u')y\\^,^Q^ Wu'h^Q) < C \\Vy (5^n")|L2(Q) lk'llL2(Q) , (4.15) 



and 



{g'^u^ - {g^u^)y) b,dy 



<C\\Vy (5V) I 



(4.16) 



with constants independent of e and ^. Applying inequahties (|4.13p - (|4.16p and then (j4.7p to (|4.12p we arrive 
at 



A, / g^v'u' dx + (A, - /3(Ae)) / g'{u'y dx 



< 



(4.17) 



< eC {\\y{gun\\h^ni) + Wa^'Whin^^uni) + ) +C^l 



■e\\2 



where C is e-independent. Since the correctors r^, rf are uniformly bounded, inequalities (|4.1ip and ()4.17p 
together imply the validity of ()3.19p . 

Step 4. Finally, it is not difficult to obtain similarly (f3?20]) via ([3lT]) . (fiTTl) and (filH]) : 



e^ao J Vu''V{g'^u^)dx 
<eC7^(||eVt;^||i2(^Q^) +e||V2^||i2(^Q^) 



< e'C l|Vn1lL2(.Qo) ||V(5'n^)|L2(,Q.,) < 



IVig^u'-' 



< 



<eCY,[h'\\LHeQo) l|V(5'^")|L2(,Q„) +£||V^/^||i2(,Qo) \\'^(.9^^')\\LHeQo)) " 



<6 (||V(5n^)||i.(^.) + < S (l|V(<7n^)||i.(^.) + ||<7n^|li2(^.uQf) + ^ 

for small enough e. 

Notice that all the estimates obtained in this section are independent of R. 



□ 



5 Some properties of two-scale convergence 

In this section we list the definitions and some properties of the two-scale convergence, see [21 [191 [211 [22] . We 
also formulate several statements (analogous to those in |21j ) which are necessary for obtaining the two-scale 
convergence of the eigenfunctions of and derivation of the limit equation. 

Let Q be an arbitrary region in R", in particular 0, = M". Denote by □ the unit cube [0, 1)". We consider 
all functions of the form u{x, y) to be 1-periodic in y in each coordinate. 

Definition 5.1. We say that bounded in Lp'iQ.) sequence is weakly two-scale convergent to a function 
V € L^{n X □), Veix) v{x,y), if 

lirn J Vi:{x)(p{x)b (^—^ dx = J J v{x,y)(p{x)b{y) dydx 
n Q □ 



13 



for all if G C^{Q) and all b G C^j.{\I\) (where C^j,(n) is the set of 1-periodic functions from C°°(M") ). 

Definition 5.2. We say that a bounded in sequence is strongly two-scale convergent to a function 

u G L?'(p, X □), ii£(x) ^ u{x,y), if 

lim J U£{x)ve:{x) dx = J J u{x,y)v{x,y) dy dx 



£- 

CI CI a 



for all Vs{x) v{x,y). 

Proposition 2. (Properties of the two-scale convergence.) 
(i) If Ue(x) u{x,y) and a G L^j.(n) then 

2 

a{x/s)ue{x) a{y)u{x,y). 

(a) Vs{x) —I- v{x,y) if and only if v six) v{x,y) and 

lim J v^dx = J J dy dx. 
ci n □ 

(itt) Iffeix) ^ f{x) m L^{n), then ^ fix). 

Proposition 3. (The mean value property of periodic functions.) Let $(?/) G Lpg^(n) . Then for each 
(/•(x) G C^(R") we have 

lim J (t)ix)^ix/e)dx = i^)y J 4>ix)dx. 

Potential vector space V^ot is defined as a closure of the set {Vtp : tp G C^^iU)} in L^(n)". We say that 
a vector b G L^(n)"' is solenoidal (6 G 14oi) if it is orthogonal to all potential vectors. Thus, 



and 

Lemma 5.3. Let and eVus be bounded in L^(M"). Then (up to a subsequence) 



u,ix)^uix,y)eL\R^,HU, 
eVueix) Vyuix,y), 

where LLp^r — -^perC-') ^^^^ Sobolev space of periodic functions. 
Lemma 5.4. Let G H'^iW'), 

Ueix)^uix)^H^iM^), (5.1) 
and «s bounded in L^(R"). Then, up to a subsequence, 

Vusix) ^ Vuix) + vix, y), where v G L^iW\ V^ot)- (5-2) 
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Lemma 5.5. Let \5.1\) and \5.0i) he valid. Let also 



for any ip e C^{Q.i) and w € C^j.(n). Then the following weak convergence of the flows takes place: 

aieQ,{e-^x)Vu,{x) ^^°™Vu(2;) in Qi, 

where homogenized matrix A^°'^ is defined by 112. 10\) . 

The proofs of the hsted statements repeat the proofs of the corresponding assertions in [21] with no or 
only small alterations, and are not given here. 



Definition 5.6. LetA^, e > 0, and Aq be non-negative self-adjoint operators inL'^(W^) andTCo C L'^(W^xQ), 
see {2. 'T\j , respectively. We say that A^ — > Aq in the sense of the strong two-scale resolvent convergence if 
{Ae + ly^ fe ^ (Ao + ly^ fo as long as fe ^ /o- 

6 Strong two-scale convergence of the eigenfunctions and multiplicity of 
the eigenvalues of 

In this section we will show that the normalized eigenfunctions are compact in the sense of strong two- 
scale convergence. Namely, provided — > Aq, a sequence of normalized eigenfunctions of the operator 
A^ strongly two-scale converges, up to a subsequence, to a function u^{x,y). This implies that vP{x,y) is 
an eigenfunction corresponding to the eigenvalue Aq of the limit operator Aq. This, together with results of 
|16j . establishes an 'asymptotic one-to-one correspondence' between isolated eigenvalues and corresponding 
eigenfunctions of the operators A;, and ^o- 

Theorem 6.1. Under the assumptions of Theorem \3.1\ An is an eigenvalue of the operator Aq. Moreover, 
there exists a subsequence e such that eigenfunctions of the operator A^ strongly two-scale converge to an 
eigenfunction u^{x,y) of Aq corresponding to the eigenvalue Xq. 

Proof. Step 1. In order to establish strong two-scale convergence of the eigenfunctions = u'^ -\- we 
establish it for each of its components separately. From (j3.22p and (j3.23p it follows that 

I|S1Ihi(R"\b«) ^ C-e-"^ (6.1) 

with C independent of e and R. From this one can easily conclude that 2^ is weakly compact in i7^(M") and 
strongly compact in L^(]R"). Indeed, since are bounded in H^{M.'^) uniformly in e, 

^ no in //^(M"), (6.2) 

up to a subsequence. For any fixed R function 2^ converges to uq weakly in H^(Bji) and, hence, strongly 
in L'^{Bji) up to a subsequence. Considering a sequence of balls Br, i? G N, one can use the method of 
extracting a diagonal subsequence such that 

Uq in L'^{Br) (6.3) 

for any R > 0. 



15 



For any 6 > we can choose R such that ||ifo||L2(]Kn\^^) < S/3 and ||'U^||i2(iRn\5^) < 5/3 for sufficiently 
small £ (the latter follows from (|6.ip ). From (|6.3p it follows that ||mo — u^\\l^{Bii) < "^/^ for sufficiently small 
e. Then, up to a subsequence, 

Iko - u'^WL^iR") < 11^0 - u'^\\l'^(Br) + I1^o||l2(R"\Bjj) + 11^^ IIl2(M"\Bjj) < ^ 

for small enough e. Hence, up to a subsequence, we have 

^ uo in l2(M"). 

Then from properties of the two-scale convergence we conclude that 

^ no. (6.4) 

Step 2. Now let us consider . We denote by ff and its restrictions Velog and Vel^^ respectively, 
extended by zero to the rest of M". 

Lemma 6.2. The following convergence properties are valid for vf (up to a subsequence): 

vfix) ^ v{x,y) e L\ni,H^{Qo)), 
eVvl{x) Vyv{x,y), 

where v{x,y) is a solution to the following problem: 

- aoAyV - Xqv = XqUo, y G Qq. (6.5) 

Here uq is a function from 

Proof. Function E H^{Qq) satisfies the following differential equation: 

- e^aoAvf - Xevl = Xeu" in (6.6) 



The right hand side of this equation is of the form XgOn^u'^. By ()6.4p and the properties of the two-scale 
convergence we have 

X,eni^{x)u'{x) ^ XoeQ,{y)enA.x)uo{x). (6.7) 
Following [21j we consider more general problem 

z,eH\nf,), -e^aoAze - Xez, = f„ f, e L^{ni). (6.8) 

(It is imphcit that = = in M^X^i^.) 
Proposition 4. Let 

f-{x)^f{x,y). (6.9) 

Then 

z%x) ^ z{x,y) e L^QuH^iQo)), 
eVz^x) Vyz{x,y), 
where function z{x,y) solves the following equation: 

- aoAyZ - Xqz = f, yeQo. (6.10) 
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Proof. One can easily derive an estimate for analogous to ()3.1ip . applying to (16. 8p a reasoning similar 
to those for the solution of equation (j3.6|) . This give us the weak two-scale convergence of z*^ and eVz^ via 
Lemma [531 The result follows by a straightforward passing to the limit in the integral identity corresponding 
to (|6.8p with appropriately chosen test function. The full proof could be found in [2T] and applies to the 
present situation with no alteration. □ 

The above proposition together with (|6.7p establishes a "weak" form of the statement of the lemma, i.e. 
weak two-scale convergence of vf. We now prove that the convergence is actually strong, following again [21]. 
Multiply (j6.6p and (j6.8p by z'^ and vf respectively and integrate by parts. The left hand sides of the resulting 
equalities are identical. So, equating the right hand sides, we obtain the following identity 



f^v'^ dx = Xf, j u'^z^ dx. 

By the definition of the strong two-scale convergence we have 

limX^ J u^z^dx = \Q j j UQ{x)z{x,y) dy dx. 

Multiplying (j6.5p and (j6.10p by z and v respectively and integrating by parts it is easy to see that 



Ao y J UQ{x)z{x,y)dydx = J J f{x,y)v{x,y)dydx. 

Thus, we have a convergence of the integrals: 

lirn J f^v\ dx = j j f{x,y)v{x,y) dy dx 
Hi Qi Qo 

for any weakly two-scale convergent sequence f^. Hence, by the definition, 

vl{x) v{x,y). 

□ 

Lemma 6.3. Sequence of functions converges to zero in the sense of strong two-scale convergence: 

2 

V2 ^ as e 0. 

Proof. Straightforward from (j3.18p and Proposition [2] (Hi). □ 
Combining ()6.4p with Lemmas 16.21 and 16.31 we arrive at 

u'^(x) n°(x, y) = no(x) + v{x, y), (6-11) 
where uq G H^R""), v e L'^i^i, H^{Qo)). 
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Step 3. Now it remains to show that v?{x,y) is an eigenfunction and Aq is the corresponding eigenvalue 
of the hmit operator Aq, i.e. that vP{x, y) satisfies (|2.1ip . In order to do that we need to choose appropriate 
test-function ij:^ and pass to the hmit in the integral identity 



e^ao / Vu^ • VV'^ dx + ai / Vu^ • V^)" dx + ooVii^ • V^/^^ dx+ 



+ 02 / Vn^ • VV'^ dx = Xe / u^V^ dx 



(6.12) 



corresponding to the original eigenvalue problem ()2.3|) ~ ()2.2p . Let us take 

'(/'^(x) = 4^o{x) + ip{x)b{e~^x), 

G Co-(M"), if G C^in^), h{y) E C7o-(Qo) 
and consider each term of (j6.12p separately. Let us expand the first term: 

e^ao / Vu^VV'^ dx = e^ao / VvfVil)^ dx+ 



+e^aQ j Wv" {Vipo + ^x)Vlp) dx + qq j eVv^LpVyb{e ^x) dx. 



(6.13) 



As S/u^ is bounded in L^-norm and iV-f/;^! < Ce~^ then the first term on the right hand side tends to zero. 
From (|3.11[) and the boundedness of VV'o + bV^p we conclude that the second term also converges to zero. 
Since by Lemma 16.21 eVt;^ converges two-scale weakly, from the definition of the weak two-scale convergence 
we obtain 

linie^ao J V u'^V ij:^ dx = j j V yv{x,y)ip{x)V yb{y) dy dx. (6.14) 

Let us show that convergence property ()5.3p holds for . To this end we substitute into (I6.12p a test 
function of the form ew{£~^x)ip{x), if € Cq^{Qi), w G C|^j.(n), cf. pT]. Then all the terms except, possibly, 

aiVn'^(x) • V yw{e^'^ x)ip{x) dx 



converge to zero. As a result, the above term also converges to zero. We then apply Lemma 15.41 for Ue 
replaced by . Since v!" coincides with on 0^, by Lemma |5 . 5 1 applied to the second term on the left hand 
side of (j6.12p with ij)'^ as in (I6.13P we obtain 



e^O 



£^0 



lim ai / Vn^ • VV-^ dx = lim oi / Vn^ • VV'o dx = I A^°'^Vuo ■ VV'o dx 



(6.15) 



For small enough e the function ^/^^ is equal to ipQ in $7q, so is bounded in Qq. Since / ao|Vti'^p dx 



is bounded uniformly in s and I^^q] — > as e ^ 0, we have 



aoVti'^VV'^ dx 



1/2 



< C I ao\Vu'\dx < C\ni\'^^^al^^ 



ao| Vn^P dx 



0. 



(6.16) 
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The function coincides with u'^ on 1^1 • Then, via (16. 2p we have convergence of the last term on the left 
hand side of (f6l^ : 

lirn 02 j Vu^ • Vip^ dx = hni 02 j Vvf ■ Vipo dx — a2 J Vu^ • VV'o dx ~ J ^'"0 ' ^V'o c?^;. (6.17) 

Thus, passing to the limit as e ^ on the left hand side of (I6.12p via (I6.14p - ()6.17p . and on the right hand 
side via (|6.1ip . we arrive at 

aQ J J V yV • LfV yb dy dx + j A^'^'^V uq • V iIjq dx + a2 j Vuq ■ Vipo dx = Xq J J (uq + v){'ipQ + ipb) dy dx. 

Since the space of functions from (j6.13p is dense in V (see ()2.8p ). the latter is equivalent to (j2.1ip . It follows 
from (16. lip . Proposition [2] (ii) and the normalization of that u^{x,y) ^ 0. Thus we have proved that Aq 
and u^{x,y) are respectively an eigenvalue and an eigenfunction of the operator Aq, completing the proof of 
the theorem. □ 

Remark 4. Theorem 16.11 combined with [13^ Theorem 2] implies the existence of eigenvalues of ^0 in the 
gaps of its essential spectrum, provided ^2 is large enough and/or 02 is small enough. 

2 

Remark 5. It is not hard to show that there holds the strong two-scale resolvent convergence Aq, see 

Definition 15.61 Namely, considering the resolvent equation 

A,w'+w' = r, 

where f^, and employing essentially the same arguments as above (cf. also |2H Theorem 5.1]), one can 

pass to the limit as e ^ in the weak form of the resolvent equation choosing appropriate test functions, cf. 
a^THh - I^Tm . to obtain that ^w^, with 



AqW^ = f. 

Further, arguing as in [21^ §4.3], cf. also proof of Lemma 16.21 above, one can show that the above weak two- 
scale convergence implies the strong one, i.e. as long as which means the strong two-scale 
resolvent convergence by the definition. The latter implies in particular the strong two-scale convergence of 
spectral projectors (-Pe(A) PoW if ^ is not an eigenvalue of Aq), see [20ll22], and has other nice properties, 
however it does not imply in its own the convergence of the spectra. The latter requires an additional 
(two-scale) compactness property to hold, which Theorem 16.11 provides. 

Remark 6. The function v{x,y) could be represented as a product of uo(a;)|^^ and Ao&(y), where b{y) solves 
(j2.13p with A = Aq. Then v{x,e~'^x) strongly two-scale converges to v{x,y) by the mean value property and 
the properties of two-scale convergence. Then 

„appr . ^ ._ / M^) + v{x, X / e) , X € , . 

^'''^> - \ uq{x), X^W\%, ^^-'^^ 



also strongly two-scale converges to vP{x,y). Hence it approximates the eigenfunction n'^(x): 

|2 

Il2(ir") 
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Now, using the result of Theorem 16. II we will discuss the multiplicity properties of the eigenvalues and 
Aq. Let us assume that the multiplicity of the eigenvalue Aq of Aq is m. Suppose that for a subsequence 
Efc — > there exist / (accounting for multiplicities) eigenvalues of A^, A^^^i < Xir^^2, ■ ■ ■ < ^ek,h such that 
Xek,i ~^ ■^0) i = 1, • • • Let ul*" be the corresponding eigenfunctions orthonormalized in L^(R"). It follows 
from Theorem 16.11 that there exists a subsequence km such that 

u.'-^ul i = l,...,l, 

where are eigenfunctions of Aq corresponding to Aq- In particular, due to the strong two-scale convergence, 
we have convergence of the inner products as a consequence of the convergence of norms: 

However (u^*"" , 'U^''™)i2(]un) = 5ij. Then u'-, i = 1, . . . ,1 are also orthonormal (in Tio), i-e. there exist at least 
I linearly independent eigenfunctions of Aq corresponding to Aq- Thus, I < m. 

The results presented in p£j remain also valid for the setting of the problem in the present paper, i.e. 
when the coefficients of the divergence form operator ^4^ are of the form (12. 4p . By Theorem 4.1 of p£], if Aq 
is an eigenvalue of the limit operator lying in a gap of its essential spectrum, then for small enough e, 
there exist eigenvalues (or at least one eigenvalue) of Agr such that 

\Xe,i-Xo\<Ce^/^i = l,...,l{e). 

Moreover, again by [16^ Thm 4.1], for any eigenfunction of corresponding to Aq the related u^'^'^'^ , see 
()6.18p . can be approximated by a linear combination of the eigenfunctions of A^ corresponding to A^^j, i = 
1, . . . , l{e). Since, by the above, {u'^'^^^ ,u'^'^^^)i2(^n-^ — > dij, as e ^ 0, z, j = 1, . . . , m, it is not hard to show 
that /(e) > m. Hence we conclude that there exist exactly m eigenvalues (counted with their multiplicities) 
of A^ such that 

|Ae,i - Aol < Ce^/^ i = l,...,m, 

where m is a multiplicity of Aq. In other words there is an "asymptotic one-to-one correspondence" between 
isolated eigenvalues and eigenfunctions of the operators A^ and Aq. 

7 Identity of the essential spectra of Aq and ^40, convergence of the spec- 
tra of in the sense of Hausdorff 

By definition, the Hausdorff convergence of spectra, o-{A^) — > a{AQ) as e — > 0, means that 

• for all A G ct^Aq) there are A^ G cr(j4e) such that A^ ^ A; 

• if Ae G (^{Ae) and A^ A, then A G (7{Aq). 

We remind that As and Aq denote the 'unperturbed' operators corresponding to As and Aq, see Section [2l 
It was shown in |22j that (t{As) o"(j4o) (the spectra of both As and Aq are purely essential). In [13j it is 
proved that the essential spectrum of a divergence form operator —V • a(x)V (where a{x) > 5 > is a scalar 
function) remains unperturbed with respect to the local perturbation of the coefficient a{x). Applying this 

assertion to the operator As and its perturbation As we conclude that ct{As) = (Tess(^e) cr{AQ). Let us 
assume that (t(Ao) = (Tess(^o)- Then cjess(^e) — > o'ess(^o)- In this case Theorem 16.11 together with the results 
of [16] imply the convergence of the discrete spectra in the gaps (adisci^e) <7disc(^o)) and, consequently, 
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we would have cr^A^) — > a^Ao). However, we cannot apply the result of |13j as it is stated to the case of the 
two-scale operators and Aq. In this section we prove the stability of the essential spectrum of with 
respect to the local perturbation of its coefficients, establishing thereby the missing part of the reasoning. 
We do this by direct means using the Weyl's criterium for the essential spectrum of an operator, see e.g. [7j. 

Theorem 7.1. The essential spectra of the operators Aq and Aq coincide. 

Proof. Step 1. First we describe the domains of Aq and Aq. According to the Friedrichs extension procedure, 
see e.g. |20] . a function u belongs to T>{Aq) if and only if n = uq{x) + v{x,y) G V and there exists h € Hq 
such that 

Bo{u,w) = {h,w)no 

for all € V, see (j2.7p - (|2.9p . li u = uq + v & 'D{Aq) then uq,v £ 'D{Aq). Due to the regularity properties of 
solutions of elliptic equations, uq G H-^^^ everywhere away from the boundary of 0.2- 

Operator acting in the Hilbert space Hq was described in [22J and is generated by a (closed) symmetric 
and bounded from below bilinear form Bq{u,w) on a dense subspace V of TIq, where TYq, V and Bq{u,w) 
are defined by ()2.7p -( f2T9]) with = and Oi = M". A function u belongs to domain T)[Aq) if and only if 
u = uq{x) + v{x, y) £ V and there exists h £ Hq such that 

Bq{u,w) = {h,w)^^ 

for all u; G V. U u = uq + v £ V{Aq) then no, u G V{Aq), uq G H'^{W). 

Let ^ be a self-adjoint operator with domain I'(A) acting in a Hilbert space H. By the Weyl's criterium, 

see e.g. [7j, condition A G crcss(^) is equivalent to the existence of a singular sequence u^^^ G T>{A), i.e. such 
that 

0<Ci< ||n('=)||H<C2, (7.1) 

u^''^ weakly in H, (7.2) 

{A - A)n('') strongly in H. (7.3) 

Step 2. Let A G (Tess(j4o) and u^^"^ — Uq \x) + v^^\x,y) be the corresponding singular sequence in 
T){Aq) C Hq. We want to construct on its basis a singular sequence for the operator Aq, i.e. in V{Aq) C Hq 
and satisfying properties (|7.ip - (|7.3p . First notice that the gradient of u^q^ is bounded in L^(M"). Indeed, 
from (12. 9p and (17. Sh we have 

||Vnf)||i,(i,„) < C^o(n(^),u(^)) = CA(nW,u(^))^^ +0(1) < C. (7.4) 
Let us define a cut-off function 

VkM^) = V I 7(kl - R) 



where rj G C^(M) is such that 

vit) 



1, t < 0, 
0, t > 1. 



Consider the following sequence, u^^^r]k,Rf, G V{Aq), where Rk is chosen large enough so that ||u^'^'^(l — 
Vk,R^)\\'}l^^ ^ \- This sequence obviously satisfies (17.ip regarding the operator Aq. 

Let us check property (j7.3p . The operator Aq acts on a function u G H'^iW^) C Hq as followflll, cf. [22]. 

Let 

-V • ^^°'"Vn(x) = fix) G l2(M«). 



^If It = uo{x) +v{x,y) then Aqu = h E fio implies —V ■ A'^°"^\/uo — {h)y and —aoAyV = h{x,y), y G Qo- 
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Then, by the definition of Aq, we have 

Aou{x) = |Qir'0Q,(y)/(x) e Wo. 

Note that 



\Aou\ 



Ho 



-1/2 



L2(R")- 



For u'^^\k,Rk derive 

Thus we arrive at 

(Io-A)(n('=)%,^J 

+ 



< 



Ho 



(fc) 



Ho 



+ \Qi\-''H2 



0(1) + -O 



(k) 



L2(IRn) 



L2(IRn) 



+ 



(7.5) 



Due to (jT.ip and ()7.4p the latter converges to as /c ^ oo. Hence (|7.3p holds regarding Aq. 

Now notice that if suppn fl = 0, then u € P(Ao) if and only if n G V^Aq); besides Aqu = Aqu. 
We hence next shift the supports of the elements of the sequence away from Q2 ensuring also that the new 
sequence is weakly convergent to maintain (j7.2p . Since suppr/^^^j. is a closed ball of radius R]. + k centered 
at the origin, the shift of x by '■= (Rk + 2k + diam(i72)) S, for every k, where ^ is an arbitrary unit vector 
from M", will do the job. Hence, for the given A we have constructed a singular sequence 



w'^''>{x,y) = u^^>{x + ik,y)r]k,R^{x + ik), 

satisfying all the properties (j7.ip - (j7.3p for the operator Aq. Namely, the translational invariance of Aq in x 
ensures that (17. ip and ()7.3p are satisfied. Finally, ()7.2|) follows from the pointwise convergence of w'^^'^ to zero 
as ^ 00 (since for any fixed x, w^^\x,y) = for large enough A;). Thus A G <Tess(^o)- 

Step 3. Suppose now that A G (7ess(^o) and u^^^ = u'^\x) + v^^\x,y) is the corresponding singular 
sequence. Let R be such that Q2 C Br. There are only two alternative possibilitie^: 

• There exists a sequence 5i ^ such that for any i G N 

\\u^''\l-eBn^MHo<S^ 

for all k. 

• There exist a constant M > and subsequences k{j) 00, ^ 00 as j — > 00 such that 



(7.6) 



Mj))i 



1-0B 



R+i(j) )\\Ho 



> M 



(7.7) 



for all j. 

Let (|7.6p take place. The sequence Vm[)''^ is bounded in L^{W), cf. From ([LS]) and 

L2(Rn) = WIWho, for all / G L2(R") c Hq, 



(7i 



l-Ho and let Si :— sup A^i- Then either Si —> giving (|7.6p or Si ^ yielding (|7.7|l . 
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it follows that 

4^"^ ^u(x) in L2(M"), 

up to a subsequence. The reasoning leading to this assertion is essentially identical to the one in (j6.iP "([6 



and is not reproduced. From (17.2p and the latter we conclude that 

v^^\x,y) —u{x) weakly in TIq. 

Hence, on one hand, we have 

^u, u^^^^^ {u,u)^^ = — j dx. 

On the other hand, 

(n,f(^-))^^^ = j j uv^^Uydx= (ueQ,{y),v^''^^^^^ -{ueQ,{y),u)^^ = -\Qo\ j dx. 

Comparing the last two formulas, conclude that at n = 0, i.e. 

^ in L2(R"). (7.9) 
Denote Aqu^^'^ by g^''\x,y) = g^''\x) + h^''\x , y) G Tio- From ([73]) and ([TS]) we get the following convergence: 

Wdo'^ - a4''''iu2(R") o> 

||/,(A^)_A^W||„^^0. (7.10) 

Then from (|7.9|) we have 

g^^^ ^OinL'^iW). (7.11) 
Analogously to [21] we define a self-adjoint operator Ay acting in L'^{fli x Qo) by 

AyV = -aoAyV = p, p e L^($7i x Qo)- 

The domain of the operator, ^{Ay) C L'^{Q,i, Hq{Qq)), is the set of all the solution of this equation. Similarly 
can be defined operator Ay acting in L^(M"' x Qo). One can easily check the following properties: 'D{Ay) C 
P(^o), T^{Ay) C ViAo), cr{Ay) = a{Ay), a{Ay) C iTess(^o) and, in particular, 

(^{Ay) C a^ssiAo). (7.12) 
It is not difficult to see (by analyzing (j2.9p . see also [21]) that 

Ayv^'^') = gf^enAx)9Q,{y) + h^'l (7.13) 

(Note that AyV^''"' / ^o^^^''^-) 

Combining (j7.10p . ()7.1ip and (j7.13p we arrive at 

\\{Ay - X)v^'^\L^^n,xQo) = ll^f ^Qo(y) + h^'^ - Xv^''^\\L^in,>cQo) - 0. 

This implies that A belongs to the spectrum of Ay (notice that (|7.ip holds for v^''^ via (|7.9p ). Hence A G 
(^o), see ([732]). 
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Now let (17. 7|) hold. Consider a sequence w^^'^ = u^''^^^\l — rji^^j-^ n) € ^{Aq) (we remind that R is large 
enough to ensure CC Br). Then 




i.e. (I7.ip is satisfied for w^^\ Since the sequence 1 — 'r]i{j),R tends to pointwise, (I7.2p is valid. Analogously 
to (|7.5p we derive 

||(1o-AM^-)!|^^ = ||(Ao-Ay^-)||wo -0, (7.14) 
yielding ()7.3p . Thus, we conclude that A G cress(^o)) completing the proof of the theorem. □ 
Remark 7. Theorem 17.11 combined with [22j implies that cJess 

(Aq) = {A : /5(A) > 0} U cr{Ay). Using the 
methods of [22] it is not hard to show further that (Tess(^o) contains no point spectrum (in particular, no 
embedded eigenvalues) except if A is an eigenvalue of Ay corresponding to an eigenfunction with zero mean. 
It is natural to conjecture (cf. [22]) that, outside these eigenvalues, the spectrum is absolutely continuous 
and the "eigenfunctions of the continuous spectrum" are u{x,y,X) = uo(x,A)(l + A6(y, A)), where uo(x. A) 
are solutions of the appropriate scattering problems: 

V • A'^^^Vuo + /3(A)no = 0, x G M"\172, 

(7.15) 

a2Auo + Xuq = 0, X € 

with the appropriate matching condition at di}2 and radiation condition at infinity. A detailed study of this 
as well as of the convergence of the related generalized eigenfunctions (cf. [22] for the defect-free case) is 
beyond the scope of the present paper. 

Summarizing the main results of the present paper we conclude that Theorems 16.11 and 17.11 together with 
the results of [131 [16] (see the discussions at the end of Section [6] and in the beginning of the present section) 
establish the validity of Theorem 12.11 
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